AEDC-TR-76-31 


•  •  .  APR  1  197? 


ffff  0  9  mi 


STATUS  REPORT  ON  THE  CALCULATION  OF  BLUNT  BOOT 
FLOW  FIELDS  IN  SUPERSONIC,  NONUNIFORM  FLOWS 


PROPULSION  WIND  TUNNEL  FACILITY 
ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
AIR  FORCE  SYSTEMS  COMMAND 
ARNOLD  AIR  FORCE  STATION,  TENNESSEE  37389 


March  1976 


Final  Report  for  Period  1  July  1974  —  30  June  1975 


Approved  for  public  release;  distribution  unlimited. 


*£^0C  LitirtAHV 
F40600- 76 -C -0001 


Prepared  for 


DIRECTORATE  OF  TECHNOLOGY 

ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 

ARNOLD  AIR  FORCE  STATION,  TENNESSEE  37389 


•  % 


NOTICES 


When  U.  S.  Government  drawings  specifications,  or  other  data  are  used  for 
any  purpose  other  than  a  definitely  related  Government  procurement 
operation,  the  Government  thereby  incurs  no  responsibility  nor  any 
obligation  whatsoever,  and  the  fact  that  the  Government  may  have 
formulated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  implication  or 
otherwise,  or  in  any  manner  licensing  the  holder  or  any  other  person  or 
corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use,  or 
sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


Qualified  users  may  obtain  copies  of  this  report  from  the  Defense 
Documentation  Center. 


References  to  named  commercial  products  in  this  report  are  not  to  be 
considered  in  any  sense  as  an  endorsement  of  the  product  by  the  United 
States  Air  Force  or  the  Government. 


APPROVAL  STATEMENT 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


FOR  THE  COMMANDER 


ELTON  R.  THOMPSON 
Research  &  Development 
Division 

Directorate  of  Technology 


ROBERT  0.  DIETZ 
Director  of  Technology 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  report  number  2  GOVT  ACCESSION  NO. 

AEDC-TR- 76-31 

3  RECIPIENT'S  CAT  ALOG  NUMBER 

4  TlTL£  fond  SubltlfeJ 

STATUS  REPORT  ON  THE  CALCULATION  OF 

BLUNT  BODY  FLOW  FIELDS  IN  SUPERSONIC, 
NONUNIFORM  FLOWS 

S  TYPE  OF  REPORT  A  PERIOD  COVERED 

Final  Report,  1  July 

1974  -  30  June  1975 

6  PERFORMING  ORB.  REPORT  NUMBER 

7  AUTHORfsJ 

Harold  J.  Schmidt,  ARO,  Inc. 

B  CONTRACT  OR  GRANT  NUMBERrO 

9  PERFORM'NG  organization  NAME  AND  ADDRESS 

Arnold  Engineering  Development  Center  (DY) 
Air  Force  Systems  Command 

Arnold  Air  Force  Station,  TN  37389 

11  controlling  office  name  AND  ADDRESS 

Arnold  Engineering  Development  Center 
(DYFS) ,  Air  Force  Systems  Command, 

Arnold  Air  Force  Station,  TN  37389 

12  REPORT  OATE 

March  1976 

13.  NUMBER  OF  PAGES 

35 

1A  monitoring  agency  name  a  ADDRESSrl/ Nirrarwif  from  Conlrolllnt  OfficaJ 

16  SECURITY  CLASS,  (ot  tbla  report) 

UNCLASSIFIED 

IS*.  DECLASSIFICATION 'DOWN  GRADING 
SCHEDULE 

N/A 

16  DISTRIBUTION  STATEMENT  ro/ fhJ*  Repe/<> 

Approved  for  public  release;  distribution  unlimited. 

17  DISTRIBUTION  STATEMENT  (et  the  abafracf  enfared  In  Block  20,  II  dlllcrmi  Irom  Report) 

IB  SUPPLEMENTARY  NOTES 

Available  in  DDC. 

19  KEY  WORDS  (Continue  on  rovoroo  aido  H  nocmaaary  and  Idantffy  by  block  number) 

nonuniform  flow  fluid  flow 

supersonic  flow  computations 

blunt  bodies 

20  ABSTRACT  fContInum  on  rmvorae  aide  H  neeaaaary  and  Identity  by  block  number) 

It  is  difficult  to  eliminate  flow  nonuniformities  in  facilities  for 
the  testing  of  ablative  heat  shield  materials.  To  Investigate  the 
potential  effects  of  free-stream  nonuniformities,  a  modification 
of  an  existing  analytical  method  developed  by  Moretti  and  co-workers 
for  uniform  blunt  body  flows  has  been  used.  The  results  of  the 
study  indicate  that  for  the  flow  model  used,  the  free-stream 
nonuniformity  can  significantly  alter  the  flow  properties  of  the 
body  surface.  In  particular,  free-stream  Mach  number  gradients 

DD  ,^Sr73  1473  EDITION  OF  1  NOV  6S  IS  OBSOLETE 


UNCLASSIFIED 


UNCLASSIFIED 


20.  ABSTRACT  (Continued) 

can  result  in  velocity  gradients  at  the  body  surface  which,  in 
turn,  will  Influence  local  heating  rates. 


UNCLASSIFIED 


AEDC-TR-76-31 


PREFACE 
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author  of  this  report  was  Harold  J.  Schmidt,  ARO,  Inc.  The  manu¬ 
script  (ARO  Control  No.  ARO-PWT-TR-75-109)  was  submitted  for 
publication  on  June  30,  1975, 


1 


AEDC-TR-76-31 


CONTENTS 


P!-.ge 


1.0  INTRODUCTION .  5 

2.0  OUTLINE  OE  THE  METHOD 

2 . 1  Field  Point  Computations  .  9 

2 . 2  Shock  Points . 15 

2.3  Body  Points . 19 

2 . 4  Outer  Boundary . 21 

3.0  RESULTS  OF  SAMPLE  CALCULATIONS . 21 

4.0  CONCLUSIONS . 31 

REFERENCES . 33 

/ 


ILLUSTRATIONS 


Figure 

1.  Flow  Geometry  and  Coordinate  System  Definition  Used  in 

the  Analysis  of  Blunt  Body  Flow  with  a  Nonuniform  Free 
Stream .  ^ 

2.  Computational  Regions  in  the  Physical  Space  and  the 

Continuous  and  Discretized  Transformed  Space  ....  8 

3.  Diagram  of  Geometry  of  the  Shock  Boundary  Calculation 

in  the  Two-Dimensional  Case  or  a  Symmetry  Plane 
Calculation  in  the  Three-Dimensional  Case . 16 

4.  Shock  Path  in  the  Quasi-One -Dimensional  Unsteady  Flow 

Plane  Showing  the  Three  Characteristic  Directions  ...  18 

5.  Body  Path  in  the  Quasi-One -Dimensional  Unsteady  Flow 
Plane  Showing  the  Three  Characteristic  Directions  with 

the  Condition  that  u  =  0  at  the  Body  Surface .  20 

6.  Three-Dimensional  Flow  Field  for  an  8-deg  Sphere 

Cone  at  3  deg  Angle  of  Attack  in  a  Mach  Number  7.  0 
Uniform  Flow . . . 33 

7.  Nonuniform  Free-Stream  Model  Used  to  Illustrate 

Effects  of  a  Radially  Increasing  Flow  Mach  Number  ...  24 


3 


AEDC-TR-76-31 


Figure  Page 


8.  Calculated  Shock  Shape  for  a  Nonuniform  Free-Stream 
Flow  in  Which  the  Mach  Number  Radially 

Increases . 25 

9.  Calculated  Static  Pressure  Distributions  for  a  Uniform 

Free  Stream  and  Nonuniform  Free  Stream  in  Which  the 
Mach  Number  Radially  Increases . 26 

10.  Calculated  Flow  Velocities  Near  the  Body  Stagnation 

Point  for  a  Uniform  Free  Stream  and  a  Nonuniform 
Free  Stream  in  Which  the  Mach  Number  Radially 
Increases . 27 

11.  Nonuniform  Free-Stream  Model  Used  to  Illustrate 
Effects  of  a  Radially  Decreasing  Flow  Mach 

Number . 28 

12.  Calculated  Shock  Shape  for  a  Nonuniform  Free-Stream 

Flow  in  Which  the  Mach  Number  Radially  Decreases  .  .  29 

13.  Calculated  Static  Pressure  Distributions  for  a  Uniform 
Free  Stream  and  a  Nonuniform  Free  Stream  in  Which 

the  Mach  Number  Radially  Decreases . 30 

14.  Calculated  Flow  Velocities  near  the  Body  Stagnation 

Point  for  a  Uniform  Free  Stream  and  a  Nonuniform 
Free  Stream  in  Which  the  Mach  Number  Radially 
Decreases . 31 

NOMENCLATURE . 35 


4 


AEDC-TR-76-31 


1.0  INTRODUCTION 


Simulation  of  conditions  in  ground  testing  of  high  performance 
reentry  vehicles  continues  to  be  a  formidable  task.  The  high  heating 
rates  and  pressure  levels  achieved  in  flight  tests  are  difficult  to 
duplicate  in  ground  tests,  where  not  only  the  test  model  but  the  test 
equipment  as  well  must  withstand  the  environment.  To  date,  the  arc 
heater  has  been  the  most  successful  device  for  producing  an  environ¬ 
ment  for  testing  ablative  materials  for  reentry  vehicles.  However, 
the  arc  heater  has  disadvantages  in  that  the  flow  Mach  number  is  low, 
the  plasma  jet  enthalpy  distribution  is  nonuniform,  and  testing  is 
limited  to  small  models.  One  method  used  to  extend  the  utility  of 
the  arc  heater  is  to  envelop  the  arc  heater  flow  with  a  uniform,  high 
Mach  number,  low  enthalpy  flow  (i.  e. ,  a  cold  shroud).  The  concept 
is  that  the  high  enthalpy  flow  from  the  arc  will  be  injected  into  the 
boundary  layer  around  the  model  and  will  correctly  simulate  the 
heating  rates  while  the  surrounding  high  Mach  number  cold  flow  will 
provide  the  correct  pressure  distribution  simulation  on  the  model. 

Whether  conventional  or  shrouded,  the  arc  heaters  in  current 
use  produce  nonuniform  though  axially  symmetric  test  jets.  To 
confidently  apply  data  from  such  a  test  to  the  design  of  considerably 
scaled-up  flight  vehicles,  one  must  ascertain  the  influence  of  the  flow 
nonuniformities.  The  object  of  the  present  effort  is  to  extend  existing 
analytical  techniques  for  blunt  body  flow-field  analysis  to  include  the 
effects  of  the  nonuniform  free  stream. 

The  most  successful  analytical  method  existing  was  developed  by 
Moretti  and  his  co-workers  (Refs.  1  through  3).  This  is  a  time- 
dependent  technique  in  which  the  blunt  body  flow-field  solution  is 
taken  to  be  the  asymptotic  limit  of  an  unsteady  flow  in  which  steady- 
state  boundary  conditions  are  imposed.  The  Moretti  technique  differs 
in  two  significant  aspects  from  other  methods  of  this  same  general 
class.  First,  the  Moretti  technique  is  direct,  with  the  body  shape 
specified  and  the  shock  shape  calculated  as  part  of  the  solution,  and 
second,  the  shock  wave  is  treated  as  a  discontinuity  located  at  the 
upstream  boundary  of  the  computational  region. 

There  are  several  versions  of  the  Moretti  program  in  use,  and 
the  versions  differ  by  the  geometry  considered,  the  main  coordinate 
systems  employed,  and  the  particular  computer  for  which  they  were 
written.  For  this  study,  two  of  these  automated  versions  were  used. 
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The  first  is  a  simple  model  applicable  to  axiS3rmmetric  blunt  bodies 
at  zero  angle  of  attack.  This  version  used  a  combination  of  two- 
dimensional  cartesian  and  polar  coordinate  systems.  Due  to  its 
relative  simplicity,  this  program  was  used  to  gain  familiarity  with 
the  Moretti  technique  and  to  perform  preliminary  calculations  with 
a  nonuniform  free  stream.  The  second  version  of  the  program  is 
for  an  axissmimetric  body  at  angle  of  attack  and  uses  a  spherical 
coordinate  system.  This  second  version,  as  obtained  from  Sandia 
Corporation,  was  written  for  a  CDC  6600  computer  and  included  real 
gas  effects.  When  this  version  was  converted  to  the  IBM  system,  the 
real  gas  capability  was  deleted  in  order  to  reduce  the  storage  require¬ 
ments  and  running  time. 

While  both  of  these  programs  have  a  capability  of  extending  the 
computations  downstream  of  the  nose  region  into  the  supersonic  flow 
field,  this  option  was  not  exercised.  This  extension  is  quite  straight¬ 
forward  but  does  require  reprogramming  to  redefine  the  variables. 

The  nose  region  of  the  blunt  body,  however,  is  of  primary  interest  in 
regard  to  the  intended  application,  and  this  effort  has  concentrated  on 
calculating  this  flow  field.  Also,  since  the  zero  angle  of  attack  is  a 
special  case  of  a  more  general  angle -of- attack  calculation,  only  the 
general  case  will  be  considered  here. 

In  Section  2.0,  the  method  is  outlined.  In  Section  3.0,  the  results 
of  the  sample  computations  are  presented.  In  Section  4.0,  the  princi¬ 
pal  conclusions,  limitations,  and  possible  extensions  of  the  analysis 
are  discussed.  In  particular,  the  results  indicate  that  flow  nonuni¬ 
formities  can  alter  the  heating  rates  to  the  extent  that  they  must  be 
taken  into  account  when  scaling  experimental  data  for  different  sized 
models. 


2.0  OUTLINE  OF  THE  METHOD 


The  flow  geometry  incorporated  in  the  blunt  body  program  with  a 
nonuniform  free  stream  is  shown  in  Fig.  1 .  For  simplicity,  the  body 
shape  is  taken  to  be  an  axisymmetric,  spherically  blunt  cone.  While 
the  program  is  written  for  axisymmetric  conical  shapes,  in  principle, 
the  body  need  only  be  planar  symmetric  in  the  pitch  plane,  and  the 
blunt  nose  does  not  have  to  be  spherical.  The  primary  spherical 
coordinate  system  is  located  on  the  body  axis  and  is  aligned  with  re¬ 
spect  to  the  wind  vector.  The  origin  of  the  coordinate  system  is 
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neither  fixed  nor  required  to  be  at  the  center  of  curvature  of  the  blunt 
nose.  However,  as  a  fundamental  assumption,  the  free-stream  stag¬ 
nation  streamline  is  taken  to  be  defined  by  the  wind  vector  passing 
through  the  coordinate  origin.  For  a  spherical  nose  at  small  angles 
of  attack,  with  the  sonic  surface  occurring  in  the  nose  region,  this  is 
a  reasonable  approximation.  In  the  spherical  coordinate  system,  the 
body  centerline  and  the  x  and  z  axes  are  in  the  same  plane.  The 
meridian  angle,  4) ,  which  has  values  between  0  and  tt  in  the  leeward 
and  windward  planes,  respectively,  is  measured  in  the  x-y  plane, 
which  is  orthogonal  to  the  x-z  plane.  The  azimuthal  angle,  0.  (mea¬ 
sured  in  the  r-z  plane),  is  normally  confined  to  the  range  w/2  <  d  <  w 
but  can  be  further  restricted  for  a  particular  problem.  Note  that 
when  the  appellations  "meridian"  and  "azimuth"  are  applied  to  the 
coordinate  angles,  it  is  the  coordinate  system  that  is  referenced  and 
not  directions  on  the  body  surface. 

The  free-stream  flow  is  taken  to  be  parallel  with  a  uniform  static 
pressure.  The  velocity  and  density  in  the  free  stream  are  taken  to 
have  axially  symmetric  profiles  with  respect  to  the  previously  defined 
stagnation  line,  which  is  an  extension  of  the  coordinate  system  base¬ 
line.  The  distribution  of  the  velocity  is  depicted  in  Fig.  1. 


Figure  1.  Flow  geometry  and  coordinate  system  definition  used  in  the  anaiysis 
of  blunt  body  flow  with  a  nonuniform  free  stream. 
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The  actual  flow  computational  region  of  interest  in  the  blunt  body- 
problem  is  the  region  between  the  bow  shock  wave,  the  shape  and 
position  of  which  are  determined  as  part  of  the  solution,  and  the  speci¬ 
fied  body  surface .  Closure  of  the  computational  region  is  pro-vided  by 
the  specification  of  a  limiting  -value  of  the  azimuth  angle,  9-may  This 
angle  is  selected  in  conjunction  with  the  placement  of.  the  coordinate 
origin  such  that  the  flow  through  this  closure  surface  is  everywhere 
supersonic.  That  is,  the  subsonic  flow  in  the  stagnation  region  of 
the  blunt  body  must  be  wholly  embedded  within  the  computational 
region.  The  computational  region  in  the  physical  space  is  depicted 
in  Fig.  2a.  The  solution  on  the  closure  surface,  or  immediately 
below  it,  will  permit  the  continuation  of  the  flow-field  calculation 
further  downstream  using  conventional  steady- state  analytical  methods 
such  as  the  method  of  characteristics. 


a.  Calculatlonal  region  in  physical  space 


b.  Computational  region  in  c.  Discretized  computational  regions 

transformed  space  in  transformed  space 

Figure  2.  Computational  regions  in  the  physical  space  and  the  continuous 


and  discretized  transformed  space. 
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The  basic  concept  is  to  consider  the  time- dependent  equations  of 
motion  for  an  inviscid  fluid  within  the  outlined  region  of  interest.  In 
a  steady-state  formulation  of  the  problem,  the  describing  equations  in 
the  subsonic  regions  of  the  flow  are  elliptic,  and  the  problem  would  be 
an  elliptic  boundary  value  problem  with  unknown  free  boundaries, 
namely  the  shock  and  sonic  surfaces.  In  the  unsteady  formulation  the 
appropriate  flow  equations  are  hyperbolic  in  character  over  the  entire 
flow  region,  and  the  problem  can  be  formulated  as  an  initial  value 
problem,  which  is  not  as  difficult  a  computational  problem.  While 
the  body  surface  is  fixed  in  time,  the  boundary  represented  by  the 
shock  surface  is  allowed  to  move,  and  in  essence  it  is  the  asymptotic 
time  position  of  the  shock  wave  which  represents  the  final  solution. 
The  calculations  in  the  interior  of  the  flow  field  are  different  from 
those  performed  at  the  boundaries,  and  each  of  these  calculations 
will  be  considered  separately  in  the  following  subsections. 


2.1  FIELD  POINT  COMPUTATIONS 


With  U0,  Up,  and  u^  representing  the  velocity  components  in 
the  spherical  coordinate  directions,  the  governing  system  of  non¬ 
steady  flow  equations  for  an  inviscid,  nonconducting,  perfect  gas  is 
given  as 


dp 

di 


a  ,  V  2  1  d  ,  . 

+  _(pup  +  -pu^  +  ^(pu^)  + 

^  da^  110  du^ 

dl  ^  dr  ^  t  dd  ^  r  sin  6  d<f> 


r  tan  0  r  dd  ^ 

“  —  r  _  _ 


0 


0 


(1) 

(2) 


dug 

dt 


+  U 


du 

^~d 


6  “0 
+  — 


du0 

dd 


"r"d 


r  sin  0  d<f> 


+  cot  0  +  —  ^ 
r  pt  dd 


=  0 


du 


“0 

+  U  -jr  ■  + - —  +  - - 

'■  dr  r  00  r  sin  0  d<f> 


0  ^ 


ds 

dt 


+ 


".^“0 


cot  0  + 


ds  "0  ds 

''•'dr  ^  r  sin  0  drf) 


1  dp^ 
pr  sin  0  d(f> 


dug 

w 


1,  "“0 
+  -  U/1 
r 


0 

0 


(3) 

(4) 

(5) 
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where  p  is  the  static  density,  p  the  static  pressure,  and  s  is  the 
entropy.  The  first  equation  is  the  continuity  equation,  followed  by  the 
three  Euler  equations,  and  the  final  equation  expresses  the  fact  that 
the  entropy  is  constant  along  the  streamlines.  In  view  of  the  complex 
geometry  involved  in  the  flow  field  and  the  fact  that  there  is  no  length 
scale  associated  with  an.inviscid  flow,  it  is  convenient  to  transform 
the  independent  variables  and  nondimensionalize  the  dependent  vari¬ 
ables  in  the  above  equations. 

To  transform  the  independent  variables,  one  can  take  the  body 
surface  as  an  explicitly  defined  function,  B,  which  is  dependent  on 
the  azimuth  angle  and  meridian  angles  only.  The  body  surface  is 
defined  by  an  equation  of  the  form 


r  -  B(0,0)  =  0 


(6) 


A  similar  relation  can  be  defined  for  the  shock  wave  surface  as 


r  -  S(^,0,l)  =  0  ^>7^ 

where  the  function  S  is  implicit  and  also  functionally  time  dependent 
since  the  shock  is  allowed  to  have  a  radial  motion.  With  5  represent¬ 
ing  the  distance  between  the  body  and  shock  surfaces  along  any  radius, 
this  distance  is  given  as 


5(<6,0.t)  =  S  -  B  (8) 

A  transformation  of  variables  can  then  be  defined  by 

c.  i^,0.  e,*  .  ^.T  .  . 

where  0,  and  T  are  the  new  independent  variables.  Under  this 
transformation  the  physical  space  of  Fig.  2a  can  be  considered  to 
map  into  the  space  of  Fig.  2b  where  the  shock  surface  boundary  lies 
in  the  plane  ?  =  1  and  the  body  surface  lies  in  the  plane  ?  =  0.  The 
required  transformation  relations  are  given  as 

d  _  I  d 
d7~  s'^ 

d_  _  ^ 

dd  ~  S  dC  ^  dQ 
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d  C($)  d  d 

^  ~  ~  S  dZ 

d  _  ids  d  ^  ± 

57  “  “  S 


(10) 


with  the  function  C(x)  being  defined  as 


C(x)  =  (1  -  o  - 

ax 

where  x  is  a  dummy  notation  representing  either  (t)  or  0  as  the  case 
may  be.  With  the  establishment  of  the  transformation  relations,  the 
original  notations  6,  0  ,  and  t  will  be  retained  for  the  transformed 
variables  to  which  they  are  identical. 


/■il 

-5x  (11) 


To  nondimensionalize  the  dependent  variables,  it  is  convenient 
to  use  the  free-stream  static  pressure  and  free-stream  density  on 
the  principal  coordinate  axis.  The  following  nondimensional  variables 
are  introduced: 


Pressure 

P  =  ?n(p/pj 

Density 

R  =  Mp'pJ 

Entropy 

E  =  P  -  yR 

Velocity 

(u,v,w)  =  (u^,  ug. 

(12) 

A  convenient  reference  length  is  taken  as  the  body  nose  radius,  and 
the  reference  time  unit  is  taken  as  r^/ 

Incorporating  the  transformation  relation,  Eq.  (10),  and  the  non - 
dimensionalizations,  Eq.  (12),  into  the  governing  system,  Eqs.  (1) 
through  (5),  yields  the  following  coupled  system  of  equations. 


dP  dP  dP 

(9u  5v  <9v 

+  “3  aj  ^  ^  ^  ^ 

dw  (9w 

■  “6  ^  +  ‘‘S 

(13) 

da  da  da 

dP 

(14) 

“0  55  5^  “2 
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— 

Bi 


—  <jw 


-BfL 

Bi 


Bv 

dv 

BP 

BP 

Bd 

+  Qi  —  4. 

^  B6 

“2  a? 

®ii  w 

+  ®12  +  ®13 

a* 

Bw 

BP 

ap 

Be 

“1  a^ 

a 

Bo  -  -T- 

^  BC 

“14  ^ 

+  ®15  +  ®16 

BE  BE  „  5E 

®o  ae  ^  bc 


(15) 

(16) 

(17) 


Si+B 

w  1 

sin  6  5^T~B 

^2  =  ^  If  -  S^Vb 


“4 


= 

KTb 

-  -y  t:(e) 

8  (8^  +  B) 


w  C(s6)~|l 
sin  9  J/ 


(18) 


6  sin  6  +  B) 

^  -y  C(<^) 

^  5  sin  9  {8^  +  B) 

y(%'  cot  9  T  2u) 

“a  =  SC  +  B 


®io 

®ii 

®12 


-  (v2  +  w2) 

5<^  +  B 


y(8C  +  B) 

-  C(a) 

yS  (S4^  +  B) 
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®13 


“15 
“16 

This  system  of  equations  is  numerically  integrated  with  respect 
to  time  in  the  interior  of  the  transformed  computational  space  shown 
in  Fig.  2b.  The  spatial  derivatives  appearing  on  the  right-hand  side 
of  these  equations  are  approximated  by  finite  differences,  and  to  do 
this,  the  continuous  space  of  Fig.  2b  is  divided  into  a  discrete  three- 
dimensional  grid  of  points  as  indicated  in  Fig.  2c.  The  space  is 
first  sectioned  by  cutting  planes  for  equally  spaced  increments  in  the 
azimuthal  angle  ^ ,  and  up  to  nine  such  planes  can  be  used  in- the  pres¬ 
ent  program.  Each  of  these  planes  is  then  divided  into  a  matrix  of 
equally  spaced  points  in  the  0  and  t  directions,  and  up  to  six  points  in  ' 
?  and  10  points  in  9  can  be  used.  This  spacing  provides  an  indication 
of  the  complexity  of  the  machine  calculations  when  one  considers  that, 
as  a  minimum,  values  for  each  of  the  five  dependent  variables  have 
to  be  stored  for  each  of  the  possible  540  grid  points.  For  an  axisym- 
metric  body  at  zero  angle  of  attack,  the  grid  equivalent  of  a  single  q 
plane  is  required,  and  the  solution  of  this  case  would  always  be  more 
detailed  in  that  more  points  could  be  stored,  than  in  the  more  complex 
case  with  the  body  at  an  angle  of  attack. 

With  values  of  the  variables  assigned  at  each  of  the  mesh  points, 
the  system  of  equations  is  numerically  integrated  with  respect  to  the 
time  to  obtain  new  values  at  the  interior  field  mesh  points.  The 
starting  values  for  the  integration  are  assumed,  but  vary  with  time 
under  the  influence  of  the  applied  boundary  conditions. 

McCormack's  differencing  scheme  (Ref.  4)  is  used  for  the  time- 
wise  numerical  integration.  This  method  is  a  two-step  method  which 
has  the  advantage  that  while  being  of  second  order  in  accuracy,  it 
does  not  require  evaluation  of  second-order  derivatives.  The  scheme 
is  as  follows.  Let  g  represent  the  value  of  any  of  the  dependent  vari¬ 
ables,  gl"  .  represent  the  value  of  the  variable  at  a  time,  t,  at  the 

-K 


2  n 

uv  —  w  cot  (7 

+  B 


Y  sin  6  -h  B) 

-  a^_ 

Sy  sin  9  {8^  +  B) 

0  -r  2u - 5> 

(04  -I-  B)  I  sin 
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spatial  position  0j,  ?jj.  In  the  first  step,  a  first  estimate  of  the 
value  of  the  variable  at  time  t  +  At  is  provided  by 


_t+At 


+  Atgi) 


(19) 


The  derivative  in  this  equation  is  given  by  the  right-hand  side  of  the 
appropriate  equation  from  Eqs.  (13)  through  (17),  with  the  spatial 
derivatives  calculated  using  backward  differences.  Thus,  for  example, 
the  difference  approximation  to  the  derivative  9  g  /  9  0  is  taken  to  be 
equal  to 

“  ®i-],j,k)/(®i  “  ^i-l) 


In  the  second  step,  a  revised  estimate  of  the  value  of  g  at  t  +  At 
is  given  by 


i+Ai 


i.j.k 


+ 


_t+At 

^i.j.k 


+ 


(20) 


The  derivative  in  this  expression  is  calculated  as  in  the  first  step  ex¬ 
cept  that  the  estimates  are  used  and  a  forward  difference  approxi- 

i»  3  > 

mation  to  the  spatial  derivatives  is  used.  The  forward  difference 
approximation  of  9  g/ 9  0  is  given  as 


^^i+l,j,k  ^i.j.k/ 


(9..I  -  9l) 


Note  that  since  the  second  step  uses  differences  of  values  calculated 
in  the  first  step,  the  overall  storage  requirement  is  effectively  doubled. 
This  is  the  cost  that  is  paid  for  retention  of  second-order  accuracy  with¬ 
out  the  calculation  of  second  derivatives.  Earlier  versions  of  the 
Moretti  analysis  used  a  second-order  Taylor  series  expansion  for  this 
time  integration  which  did  require  second- order  derivative  calculations. 

In  performing  the  integration  with  respect  to  time,  one  cannot 
select  the  time  increment  arbitrarily,  since  it  must  satisfy  the 
Courant-Friedrichs-Lewy  (Ref.  4)  stability  criteria  given  as 
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At 


< 


( 


Ar 

u  +  a 


rA0 

u  +  a 


r  sin  0  A(^ 
u  +  a  , 


(21) 


where  Ar,  rA  ^ ,  r  sin  0  A<f>  are  the  physical  space  mesh  sizes  and  u 
and  a  are  the  total  velocity  and  acoustic  speed,  respectively,  at  a 
grid  point.  Simply  stated,  for  the  calculation  to  be  stable  the  time 
increment  must  be  selected  to  be  smaller  than  the  smallest  physical 
mesh  dimension  divided  by  the  largest  signal  propagation  speed  at 
any  of  the  grid  points.  In  the  present  method,  the  time  increment  is 
taken  to  be  two -thirds  of  the  theoretical  maximum  value. 


The  application  of  this  stability  criterion  leads  to  an  important 
observation.  If  one  should  attempt  to  obtain  a  more  detailed  flow- 
field  definition  by  halving  the  spatial  increments,  not  only  would  the 
computational  time  increase  due  to  the  eightfold  increase  in  the  num¬ 
ber  of  grid  points  (fourfold  for  the  2-D  case),  but  also  twice  as  many 
timewise  integrations  would  be  required  to  achieve  the  same  asymp¬ 
totic  solution.  It  is  also  observed  that  the  use  of  grid  spacings  which 
are  nearly  equal  is  efficient  since  the  smallest  grid  spacing  is  the  one 
that  controls  the  time  increment. 


2.2  SHOCK  POINTS 

The  points  on  the  boundaries  of  the  computational  region  are 
computed  differently,  and  the  shock  boundary  layer  through  which 
the  free -stream  nonuniformity  enters  will  be  considered  first. 
Unfortunately,  it  is  not  possible  to  illustrate  the  four-dimensional 
space  in  which  the  calculation  is  carried  out,  so  for  purposes  of  illus¬ 
tration,  the  two-dimensional  analog  will  be  used.  This  model  is  also 
applicable  to  the  angle -of- attack  computation  in  the  leeward  and  wind¬ 
ward  meridian  planes  where  the  shock  curvature,  in  directions 
normal  to  these  planes,  is  symmetric. 

To  illustrate  the  shock  calculation,  the  physical  system  is  illus¬ 
trated  in  Fig.  3  for  the  leeward  plane,  (j)  =  0.  Consider  that  at 
time  t  the  solution  in  the  physical  plane  is  known  and,  in  particular, 
a  point  on  the  shock  wave  Aq  has  a  radial  velocity,  W.  In  the  time 
interval  At,  the  shock  point  moves  a  distance  WAt  and  can  be  located 
at  the  point  A  in  the  t  +  At  plane.  The  position  of  A  relative  to  the 
origin  of  the  coordinate  system  determines  the  free -stream  flow  con¬ 
ditions  upstream  of  the  shock. 
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Line  AQ  Is  backward- running  characteristic. 

Figure  3.  Diagram  of  geometry  of  the  shock  boundary  calculation  in  the 
two-dimensional  case  or  a  symmetry  plane  calculation  in 
the  three-dimensional  case. 


At  the  point  A  let  N  be  an  inward-directed  vector  normal  to  the 
sh(^k  surface  and  define  an  auxiliary_cartesian  coordinate  system 
J  t)  at  A  such  that  the  unit  vector  e  is  parallel  to  N  and  n  is  normal 
to  ?  in  the  plane  (A  =  constant.  The  unit  vector  t  is  normal  to  €  and  rj. 
Furth^,  let  the  velocity  components  in  the  directions  T,  i? ,  t  be  given 
as  u,  V,  and  w  ,  respectively.  The  direction  t  an'd  velocity  w  are 
omitted  from  Fig,  3  for  clarity. 

It  is  well  known  that  flow  velocity  components  tangent  to  an 
oblique  shock  are  unaltered  by  the  shock  while  the  normal  component 
is  changed  in  accordance  with  the  Rankine-Hugoniot  equations. 

Defining  as  the  component  of  the  shock  velocity  in  the  direction 
N  and  letting  the  subscripts  1  and  2  represent  the  conditions  fore  and 
aft  of  the  shock  at  point  A,  the  Rankine-Hugoniot  conditions  can  be 
written  as 


''1 

“2 


Vj  ,  Wj  =  Wj 

(y  -  l)('uj  -  +  2  aj 

(y  +  l)(uj  - 
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2y(U]  -  “  (y  -  n 

rn 

ly  +  n(p2/Pi)  +  (y  — 

(y  +  1)  +  (y-  l)P2/Pi 

This  is  a  system  of  five  equations  in  the  six  unknowns;  that  is,  the 
shock  position  is  known,  but  the  shock  velocity  is  unknown.  To  pro¬ 
vide  the  additional  functional  relation  to  complete  the  set,  a  quasi- 
one- dimensional  analysis  is  employed.  Consider  a  plane  defined  by 
the  ?-t  coordinates  as  indicated  in  Fig.  3.  In  this  plane  the  flow  is 
considered  to  be  quasi-one-dimensional,  nonsteady  flow  under  the 
assumption  that  the  components  of  variable  gradients  tangent  to  the 
shock  wave  are  constant! 

The  situation  in  the  ?,t  plane  is  shown  in  Fig!  4.  When  the 
quasi-one-dimensional  analysis  is  applied,  it  is  found  that  there  are 
three  characteristics  passing  through  the  point  A,  given  by 


Pi 

Pi 


dt 


u 


(23) 


and  indicated  in  Fig.  4.  The  appropriate  characteristic  is  the  first 
one  given,  and  the  compatibility  equation  along  this  characteristic  is 
given  as 


du  a  dP^ 
dt  y  dt 


3u\  /v  5P  v*  5p 


(24) 


This  compatibility  equation  provides  the  additional  relation  which 
when  added  to  the  Rankine-Hugoniot  relations  results  in  a  tractable 
set  of  nonlinear,  algebraic  equations  which  are  solved  by  an  iterative 
procedure.  The  iterative  procedure  is  as  follows: 


1 .  Assume  a  value  of  the  shock  velocity  Wn  at  point  A 
in  the  t  +  At  solution  plane. 

2.  Using  the  assumed  shock  velocity,  evaluate  the  flow 
conditions  behind  the  shock  from  the  Rankine-Hugoniot 
relation,  and,  in  particular,  determine  a  value  for  U2 . 

3 .  Construct  the  characteristic  given  by  d? /dt  =  u  -  a  and 
trace  it  back  in  time  to  the  solution  at  time  tQ.  Let  the 
intercept  be  designated  at  Q  as  shown  in  Figs.  3  and  4. 
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4.  Interpolate  in  the  time  =  to  solution  to  determine  the 
flow  qualities  and  their  derivatives  at  Q. 

5.  Using  the  values  at  point  A  from  Step  2  and  the 
values  at  Q  from  Step  4,  integrate  Eq.  {6)  from  Q 
to  A.  This  results  in  a  new  value  of  U2  at  A, 

6.  Repeat  Steps  1  through  5  until  a  shock  velocity  is 
found  which  causes  the  calculated  values  of  u2  from 
Steps  2  and  5  to  agree. 


Figure  4.  Shock  path  in  the  quasi-one-dimensional  unsteady  flow  plane 
showing  the  three  characteristic  directions. 


Since  the  free -stream  nonuniformity  enters  into  the  solution 
through  the  shock -boundary  calculations,  some  limitations  on  the 
allowable  nonuniformities  must  be  accepted.  Obviously,  the  free- 
stream  properties  must  be  continuous.  For  example,  the  procedure 
would  be  incapable  of  handling  a  discontinuity  in  free -stream  velocity 
such  as  a  shear  layer  whose  thickness  is  much  smaller  than  the 
physical  mesh  size.  In  such  a  case  one  moving  shock  point  could 
alternately  pass  in  and  out  of  the  shear  layer  with  a  resulting  cyclical 
instability.  As  indicated  in  Section  2. 1,  a  reduction  of  the  mesh  size 
to  account  for  such  a  locally  steep  gradient  in  the  free  stream  would 
detract  from  the  convergence  due  to  the  accompanying  reduction  in 
the  allowable  time  increment. 
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It  is  assumed  that  by  using  the  quasi-one-dimensional  analysis  in 
the  shock  calculation,  one  obtains  gradients  tangent  to  the  shock  which 
are  small  in  comparison  with  the  discontinuous  changes  normal  to  the 
shock.  In  the  case  of  a  uniform  free  stream,  the  gradients  behind  and 
tangent  to  the  shock  are  due  solely  to  the  fact  that  the  shock  has  a 
slowly  varying  radius  of  curvature.  The  presence  of  free -stream  non- 
uniformities  can  change  the  curvature  of  the  shock,  and  poor  conver¬ 
gence  might  occur  if  the  free-stream  nonuniformity  is  concentrated 
where  the  shock  experiences  a  rapid  change  in  curvature  (e.  g.  ,  after 
the  intercept  with  the  sonic  line).  To  summarize,  the  free-stream 
nonuniformities  must  generally  be  continuous  and  without  large  gradi¬ 
ents  so  that  significant  variations  do  not  occur  over  lengths  compar¬ 
able  to  the  mesh  size.  The  free-stream  gradients  should  not  be  sub¬ 
stantially  larger  than  those  resulting  from  shock  curvature  in  the  case 
of  a  uniform  flow. 


2.3  BODY  POINTS 

The  calculation  of  the  body  boundary  points  is  similar  to  the  shock 
calculation,  though  simpler,  since  the  entropy  at  the  body  is  known  and 
the  normal  component  of  the  velocity  vanishes.  As  in  the  shock  compu¬ 
tation,  an  auxiliary  cartesian  coordinate  system  defined  by  the  unit 
vectors  T,  r?  ,  and  "r  is  defined  but  with  5  aligned  normal  to  the  body 
surface.  The  associated  velocity  components  in  the  coordinate  direc¬ 
tions  are  again  taken  to  be  u,  v,  and  w,  respectively.  When  the  quasi- 
one-dimensional  analysis  is  applied  at  the  body,  the  appropriate  char¬ 
acteristic  is  the  one  defined  by 


in  the  ?-t  plane,  as  illustrated  in  Fig.  5.  While  u  is  known  to  vanish 
at  A,  the  acoustic  velocity,  a,  is  not  known  and  must  be  initially 
guessed  to  construct  the  backward- running  characteristic.  Assuming 
a  value  for  the  acoustic  velocity,  the  characteristic  is  constructed 
back  to  the  previous  time  solution  at  the  point  Q  where  the  flow  vari¬ 
ables  are  interpolated  and  transformed  to  the  ?,  tj  ,  t  coordinate  sys¬ 
tem.  The  appropriate  compatibility  equation  along  the  u  +  a  character¬ 
istic  is  given  as 


du  a  dP  /_  5u  _  /  v  dP  w  (9P  dv  (?w\ 

^  y  dT  \  ^  ^  dr }~  \y  d-q  y  dr  dq  ^  (26) 
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The  compatibility  equation  is  approximated  by  a  finite  difference  equar 
tion  along  the  characteristic  with  the  right-hand  side  evaluated  at  the 
point  Q.  This  evaluation  provides  a  value  of  the  pressure  P  at  point  A. 
The  calculated  pressure  and  the  known  entropy  determine  the  density 
at  A  as  well  as  a  new  value  for  the  acoustic  velocity.  An  iterative 
procedure  is  performed  until  the  resulting  calculated  acoustic  velocity 
agrees  with  the  value  used  to  construct  the  characteristic. 


d| 


a?  ■ 


0 


Figure  5.  Body  path  in  the  quasi-one-dimei»ional  unsteady  flow  plane 
showing  the  three  characteristic  directions  with  the  condition 
that  u  =  0  at  the  body  surface. 

After  convergence  of  the  iteration  to  determine  the  acoustic 
velocity  at  the  body  point,  it  is  necessary  to  determine  the  velocity 
components  v  and  w.  These  components  are  obtained  from  the  two 
equations 


and 


Jv 

di 


+  5  .  0 

P  dr) 


(27) 


(9w  _  (9w  _  ffw  p  dP  _ 

-  =  '•'  '=r~  +w=+-^=0 

ot  dT]  dr  p  ovi 


(28) 


again  approximating  the  equations  by  finite  differences.  Equations 
(27)  and  (28)  are  simply  the  inviscid  flow  equations  (3)  and  (4) 
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written  in  the  rj ,  r  cartesian  coordinate  system  with  the  u  terms 
deleted  on  account  of  the  fact  that  they  vanish  at  the  surface. 


2.4  OUTER  BOUNDARY 

One  does  not  have  a  priori  knowledge  of  any  of  the  boundary  con¬ 
ditions  along  the  surface  9  =  ^max  indicated  in  Fig.  2a.  As  a  conse¬ 
quence,  no  boundary  conditions  are  applied  here.  However,  since 
this  surface  is  totally  comprised  of  supersonic  flow,  only  those  dis¬ 
turbances  originating  in  the  interior  of  the  computational  region  can 
be  propagated  across  this  boundary.  Therefore,  the  points  on  this 
surface  are  calculated  by  simple  extrapolation  with  respect  to  the 
spatial  variables.  The  extrapolation  is  carried  out  for  each  time 
solution  generated,  but  the  temporal  derivatives  of  the  variables  at 
this  surface  never  actually  enter  into  the  computational  procedure. 

In  addition  to  the  outlined  boundary  value  calculations,  there  are 
certain  extraneous  conditions  which  should  be  noted.  The  computa¬ 
tional  planes  0=0  and  0  =  tt  are  planes  of  symmetry  so  that  deriva¬ 
tives  with  respect  to  d  vanish  identically  on  these  surfaces.  The 
stagnation  line  is  common  to  all  of  the  0  planes  and  need  only  be  calcu¬ 
lated  once.  The  stagnation  line  is  taken  to  be  a  line  of  symmetry  for  the 
evaluation  of  derivatives  normal  to  this  line. 

The  present  program  is  written  for  axisymmetric  bodies  gene¬ 
rated  from  conic  sections  (circle,  ellipse,  paraboloid)  with  a  conical 
aftbody.  The  extension  to  other  shapes  of  revolution  is  straightfor¬ 
ward;  however,  the  surfaces  must  have  continuous  second-order 
partial  derivatives.  The  restriction  to  axisymmetric  bodies  simpli¬ 
fies  the  calculations  at  the  body  boundary  since  the  normal  vector  to 
the  body  always  lies  in  the  meridian  planes.  As  pointed  out,  this 
occurs  only  in  the  windward  and  leeward  planes  in  the  case  of  the 
shock  computation. 


3.0  RESULTS  OF  SAMPLE  CALCULATIONS 


A  solution  to  the  blunt  body  problem  in  three  spatial  dimensions 
results  in  an  extreme  profileration  of  numerical  information  which  in 
itself  represents  a  challenge  to  condense.  It  is  possible  to  provide 
detail  maps  and  to  trace  streamlines  and  surfaces  of  constant  Mach 
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number,  pressure,  and  density,  but  these  are  not  of  particular 
interest.  However,  certain  gross  results  are  of  primary  interest; 
these  include  shock  shape,  sonic  line  location,  pressure  distribution, 
and  the  rate  of  convergence  to  the  solution.  Only  these  gross  char¬ 
acteristics  will  be  discussed  for  a  sample  of  cases  run  with  the  pro¬ 
gram. 

All  of  the  results  presented  are  for  a  spherically  blunt,  8-deg 
half -angle  cone  at  a  3 -deg  angle  of  attack  with  a  stagnation  line, 
free- stream  Mach  number  of  7.  The  results  were  obtained  for 
seven  <j)  planes  with  the  grid  in  each  ^  plane  having  a  total  of  42  points, 
including  the  boundary  points.  The  maximum  azimuthail  angle  is 
90  deg,  and  the  coordinate  system  origin  is  located  one-half  nose 
radius  behind  the  center  of  curvature  of  the  spherical  cap.  For  con¬ 
venience,  the  results  are  presented  in  a  conventional  manner  with 
the  body  azimuthal  angle  measured  from  the  stagnation  point. 

Figure  6  depicts  the  flow  field  computation  results  in  the  pitch 
plane  for  a  uniform  free  stream  and  illustrates  the  shock  shape 
evolution  with  respect  to  the  number  of  time  iterations  of  the  describ¬ 
ing  equations.  As  the  figure  illustrates,  the  initial  guess  of  the  shock 
shape  was  not  very  close;  however,  the  estimate  of  the  shock  standoff 
distance  was  quite  close.  The  reason  for  this  apparently  good  guess 
is  that  it  was  made  on  the  basis  of  data  correlated  for  blimt  cones  at 
zero  angle  of  attack  and  does  not  change  much  for  small  angles  of 
attack. 

The  results  of  Fig.  6  indicate  that  the  evolution  of  the  shock  shape 
is  quite  rapid,  with  an  indistinguishable  difference  between  the  results 
of  200  and  300  time  steps  for  the  scale  used  in  the  figure.  The  final 
shock  shape  in  Fig.  6  is  not  symmetric  with  respect  to  the  body  and 
is  closer  to  the  body  in  the  lower,  windward  half  plane. 

As  an  example  of  a  nonuniform  flow  calculation,  the  free-stream 
nonuniformity  was  taken  to  be  a  flow  having  a  centerline  Mach  number 
of  7  and  an  outer  flow  Mach  number  of  14.  The  velocity  distribution 
was  assumed  to  be  a  sine  function  over  a  distance  of  two  nose  radii 
normal  to  the  stagnation  line.  Profiles  of  the  density,  velocity,  and 
flow  Mach  number,  referenced  to  the  centerline  values,  for  the  non- 
uniform  free  stream  are  indicated  in  Fig.  7.  As  indicated  in  Fig.  7, 
the  resulting  velocity  variation  is  only  on  the  order  of  4  percent. 
However,  the  gradients  in  Mach  number  and  density  are  substantial. 

In  addition,  since  the  free-stream  static  pressure  is  constant,  the 
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indicated  variation  in  Mach  number  implies  an  order  of  magnitude 
variation  in  the  free-stream  total  pressure. 


Figure  6.  Three-dimensional  flow  fieid  for  an  8-deg  sphere  cone  at  3  deg 
angie  of  attack  in  a  Mach  number  7.0  uniform  flow. 
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Figure  7.  Nonuniform  free-stream  model  used  to  illustrate  effects  of  a 
radially  increasing  flow  Mach  number. 

In  Fig.  8,  the  resulting  shock  geometry  calculated  for  the  free- 
stream  nonuniformity  of  Fig.  7  is  presented.  The  shock  shape  of 
the  uniform  flow  calculation  (Fig.  6)  is  also  shown  in  Fig.  8  for  com¬ 
parison.  In  comparing  the  shapes,  one  finds  that  the  result  for  the 
nonuniform  flow  has  a  smaller  radius  of  curvature,  as  would  be 
anticipated  by  virtue  of  the  higher  free-stream  Mach  number.  By 
the  same  reasoning,  one  would  expect  that  the  shock  standoff  dis¬ 
tance  would  also  be  decreased  in  the  case  of  the  nonuniform  flow, 
although  the  results  indicate  the  contrary.  It  is  also  possible  to 
make  a  qualitative  comparison  of  the  sonic  lines  in  Fig.  8  with  those 
of  Fig.  6.  For  the  nonuniform  flow  model  used,  the  sonic  lines  are 
straighter.  The  sonic  line  intercepts  with  the  body  occur  at  nearly 
the  same  locations  for  both  the  uniform  and  nonuniform  cases. 
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However,  the  intercept  of  the  sonic  line  with  the  shock  occurs  closer 
to  the  coordinate  axis  in  the  case  of  the  nonuniform  flow.  This  can  be 
a  Mach  number  effect  indicative  of  a  steeper  velocity  gradient  behind 
the  shock. 


Figure  8.  Calculated  shock  shape  for  a  nonuniform  frea-stream  flow  in 
which  the  Mach  number  radially  increases. 

In  Fig.  9,  the  resulting  static  pressure  distributions  along  the 
body  surface  are  presented  for  the  uniform  and  nonuniform  compu¬ 
tations.  The  pressure  distributions  are  given  for  the  windward 
(<^  =180  deg),  side  (0  =  90  deg),  and  leeward  {<f>  =  0)  half  planes. 


25 


AEDC-TR-76-31 


Individually,  the  pressure  distributions  for  the  uniform  flow  case 
appear  reasonable  and  self-consistent  in  that  the  pressures  decrease 
monatonically  with  the  azimuthal  angle  and  result  in  higher  values  in 
the  windward  plane  than  in  the  leeward  plane.  The  increased  pres¬ 
sures  for  the  nonuniform  case  as  well  are  consistent  with  intuition 
since  the  free-stream  total  pressure  of  the  nonuniform  case  is 
larger  than  that  of  the  uniform  case. 


Figure  9.  Calculated  static  pressure  distributions  for  a  uniform  free 
stream  and  nonuniform  free  stream  in  which  the  Mach 
number  radially  increases. 

Of  particular  concern,  however,  is  the  pressure  overshoot  near 
the  stagnation  point  for  the  nonuniform  flow.  It  would  appear  that 
such  an  overshoot  is  not  valid  and  may  be  attributed  to  a  poor  definition 
of  the  flow  field  or  a  violation  of  the  assumption  that  the  stagnation  line 


26 


AEDC-TR-76-31 


lies  on  the  coordinate  axis.  On  the  other  hand,  in  support  of  the  exist¬ 
ence  of  such  an  overshoot,  the  results  for  all  the  variables  are  con¬ 
sistent.  The  velocity  distributions  in  the  vicinity  of  the  stagnation 
point  are  shown  in  Fig,  10,  and  the  indicated  reduction  in  the  velocity 
in  the  case  of  the  nonuniform  flow  appears  consistent  with  the  pressure 
increase  shown  in  Fig.  9.  In  contrast  to  the  earlier  discussion  of  the 
sonic  line  shape  for  the  nonuniform  flow,  the  results  indicate  that  the 
velocity  gradients  behind  the  shock  are  increased,  but  are  reduced  at 
the  body  as  a  consequence  of  the  nonuniformity  in  the  free  stream. 


Figure  10.  Calculated  flow  velocities  near  the  body  stagnation  point 
for  a  uniform  free  stream  and  a  nonuniform  free  stream  in 
which  the  Mach  number  radially  increases. 

If  this  is  indeed  the  case,  it  can  have  a  significant  impact  on  the 
interpretation  of  experimental  stagnation  point  heating  rates.  Stagna¬ 
tion  point  heating  rates  are  directly  proportional  to  the  square  root  of 
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the  velocity  gradient,  and  the  velocity  gradient  for  a  uniform  flow  is 
inversely  proportional  to  the  nose  radius.  This  is  not  true  for  a 
nonuniform  flow,  and  simple  correlation  of  stagnation  point  heating 
with  nose  radius  may  not  be  valid. 

The  flow  nonuniformity  can  also  be  reversed  in  the  sense  that 
the  outer  flow  Mach  number  can  be  less  than  the  centerline  value. 
Results  of  such  a  computation  with  the  outer  Mach  number  assumed 
to  be  5  are  shown  in  Figs.  11  through  14.  The  velocity  distribution 
was  assumed  to  be  a  sine  function,  as  described  previously.  In 
addition  to  the  noticeable  changes  in  the  sonic  lines,  the  shock  wave 
shape  that  results  exhibits  a  noticeable  change  in  curvature.  This 
results  from  the  nonmonatonic  variation  of  the  shock  curvature.  It 


. . . L_l 
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Velocity  Ratio,  u/u^ 


Figure  11.  Nonuniform  free-stream  model  used  to  illustrate  effects  of 
a  radially  decreasing  flow  Mach  number. 
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would  be  expected  that  if  the  outer  flow  Mach  number  were  reduced 
too  much  the  sonic  flow  wovild  extend  outward  to  the  point  where  it 
would  not  be  possible  to  retain  a  reasonable  mesh  size  and  satisfy 
the  criterion  that  the  flow  across  the  closure  surface  be  supersonic. 
Finally,  it  is  noted  that  the  solution  is,  in  fact,  opposite  to  that  of 
the  case  previously  presented  in  that  the  stagnation  point  pressures 
were  lowered  (Fig.  13)  and  the  velocity  gradients  increased  over 
those  of  the  uniform  flow  case  (Fig.  14). 


Figure  12.  Calculated  shock  shape  for  a  nonuniform  free-stream  flow 
in  which  the  Mach  number  radially  decreases. 
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Figure  14.  Calculated  flow  velocities  near  the  body  stagnation  point 
for  a  uniform  free  stream  and  a  nonuniform  free  stream  in 
which  the  Mach  number  radially  decreases. 


4.0  CONCLUSIONS 


The  primary  conclusion  to  be  drawn  from  this  study  is  that  a 
free- stream  nonuniformity  of  the  type  considered  can  indeed  have  a 
significant  effect  on  the  results  obtained  from  model  tests  in  such 
an  environment.  In  particular,  the  free -stream  effects  must  be 
properly  accounted  for  when  one  is  correlating  or  comparing  experi¬ 
mental  data  from  different  sources.  A  second  conclusion  which  the 
results  support  is  that  the  concept  of  simulating  the  pressure  distri¬ 
bution  around  a  model  by  the  use  of  a  shrouded  flow  is  viable.  The 
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results  support  the  observation  that  Moretti's  time-dependent  method, 
originally  created  as  a  design  tool  in  uniform  flows,  has  a  potential 
for  significant  utilization  in  experimental  investigations  where  the 
desirable  flow  uniformity  is  difficult  to  achieve. 

In  its  present  form,  the  programmed  analysis  restricts  the  form 
of  the  allowable  nonuniform  free -stream  definition  to  be  axially  sym¬ 
metric  about  a  body  stagnation  line.  In  addition  to  this  restriction, 
there  are  additional  restrictions  for  which  no  firm  guidelines  can  be 
established.  Thus,  one  has  to  take  into  consideration  the  interrela¬ 
tions  between  the  free-stream  Mach  number,  the  grid  dimensions, 
the  coordinate  region,  and  the  gradient  in  the  free  stream  when 
formulating  a  problem.  Many  of  these  aspects  have  been  considered 
in  the  appropriate  preceding  sections;  however,  fuller  understanding 
of  and  familiarity  with  the  constraints  comes  only  with  application  to 
problems  of  interest. 

The  flow  Mach  number  is  of  primary  importance  in  the  calculation. 
Indeed,  as  pointed  out,  the  mathematical  character  of  the  describing 
steady-state  flow  equations  can  be  elliptic  or  hyperbolic  depending  on 
whether  the  Mach  number  is  less  than  or  greater  than  unity.  Generally, 
the  numerics  of  the  computation  become  more  inherently  stable  with 
increasing  free-stream  Mach  number.  This  results  from  the  stability 
criteria  as  well  as  the  fact  that  the  flow  tends  to  be  asymptotic  in  Mach 
number.  In  other  words,  two  flow  geometries  for  free-stream  Mach 
numbers  of  10  and  12,  for  example,  are  physically  quite  similar;  more 
so  than  two  flows  having  Mach  numbers  of  2  and  3. 

For  a  uniform  free  stream,  experience  has  shown  that  the  solution 
starts  to  exhibit  poor  convergence  to  the  asymptotic  limit  as  the  Mach 
number  approaches  1.  5.  For  a  nonuniform  free  stream,  this  limit 
can  serve  as  a  guideline,  but  no  further  firm  criteria  can  be  estab¬ 
lished.  For  example,  it  can  be  expected  that  for  a  free  stream  in 
which  the  Mach  number  varies  over  the  range  from  2.0  to  2.  5,  a 
solution  can  be  readily  obtained.  On  the  other  hand,  if  the  variation 
is  over  the  range  from  2.0  to  10.0,  the  method  would  break  down  as 
a  consequence  of  the  other  large  gradients  in  the  free  stream. 

Certainly  in  the  case  of  zero  angle  of  attack  (two  dimensions)  and 
possibly  in  the  case  of  a  finite  angle  of  attack,  the  calculation  can  be 
extended  further  downstream  on  the  body  using  conventional  steady 
flow  analytical  methods.  At  the  present  time,  such  an  extension, 
while  of  academic  interest,  would  not  be  applicable  to  any  known 
practical  application. 
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The  present  method  is  also  incapable  of  accounting  for  disconti¬ 
nuities  such  as  thin  shear  layers  or  impinging  shocks  in  the  free 
stream.  Such  surfaces  would  actually  be  additional  boundary  value 
surfaces.  It  would  appear  that  such  phenomena  could  be  considered 
on  an  individual  basis,  but  such  consideration  would  require  a  com¬ 
plete  reformulation  of  the  analysis. 

A  basic  assumption  in  the  method  is  that  the  stagnation  streamline 
is  straight,  parallel  to  the  free-stream  wind  direction,  and  normal  to 
the  bow  shock  wave.  This  assumption  permits  the  steady  normal  shock 
entropy  to  be  assigned  to  the  body  surface.  In  a  uniform  free  stream, 
this  assumption  is  valid  only  if  the  body  surface  in  the  subsonic  region 
of  the  flow  field  is  symmetric  about  the  stagnation  streamline.  This 
does  not  require  that  the  body  be  axisymmetric;  however,  spherically 
blunted  shapes  which  are  of  common  interest  always  satisfy  this  require¬ 
ment  as  long  as  the  angle  of  attack  is  such  that  the  subsonic  flow  field 
does  not  extend  beyond  the  spherical  cap.  In  the  case  of  a  nonuniform 
free  stream,  the  validity  of  the  above  assumption  additionally  requires 
that  the  free-stream  properties  have  a  similar  symmetry  about  the 
stagnation  streamline.  In  the  present  calculation,  the  axisymmetric 
model  of  the  nonuniform  free  stream  is  representative  of  the  situation 
one  would  expect  to  exist  in  conventional  aerodynamic  testing  devices. 
For  more  general  tip  shapes  or  flow  nonuniformities,  the  assignment 
of  the  body  entropy  at  the  body  surface  would  be  an  approximation,  and 
results  based  on  the  method  would  have  to  be  considered  accordingly. 


REFERENCES 


1.  Moretti,  G.  and  Abbett,  M.  "A  Fast,  Direct,  and  Accurate 

Technique  for  the  Blunt  Body  Problem,  Part  I  -  Analysis. " 
Technical  Report  583,  General  Applied  Science  Laboratories, 
Inc. ,  June  1966. 

2.  Moretti,  G.  and  Bleich,  G.  "Three-Dimensional  Flow  Around 

Blunt  Bodies. "  AIAA  Journal,  Vol.  5,  No,  9,  September 
1967,  pp.  1557-1562. 

3.  Abbett,  M.J.  and  Fort,  R.  "Three-Dimensional  Inviscid  Flow 

About  Supersonic  Blunt  Cones  at  Angle  of  Attack,  III: 

Coupled  Subsonic  and  Supersonic  Programs  for  Inviscid 
Three-Dimensional  Flow.  "  Sandia  Laboratories  Report 
SC-CR-68-3728,  September  1968. 


33 


AEDC-TR-76-31 


4.  MacCormack,  R.  W.  "The  Effect  of  Viscosity  in  Hypervelocity 

Impact  Cratering.  "  AIAA  Paper  No.  69-354,  April  1969. 

5.  Courant,  R. ,  Friedrichs,  K.  O. ,  and  Lewy,  H,  "Uber  die 

partiellen  Differenzengleichungen  der  Mathematischen 
Physik. "  Math.  Annual,  Vol.  100,  p.  32,  1928. 


NOMENCLATURE 


a 

ao 
B 
C 
E 

g 

i,  j,k 
M 
P 
P 

Pco 

R 

r 

S 

s 

T 

t 

too 

U,  V,  w 
U,  V,  w 


Speed  of  sound,  cm/ sec 

aj  5  Coefficients  in  governing  flow  equations 

Function  defining  body  surface 
Auxiliary  function  defined  by  Eq.  (7) 

Nondimen sional  entropy,  E  =  P  -  TR 
Dummy  notation  for  a  dependent  variable 
Indices  for  spatial  grid  in  0 ,  0,  t  space 
Mach  number 

Nondimensional  pressure  =  J2n(p/p„(g^) 

Static  pressure,  dynes/cm^ 

Free -stream  static  pressure,  dynes/ cm'^ 
Nondimensional  density  =  -^nlp/p^^) 

Radial  coordinate,  cm 

Body  nose  radius,  cm 

Implicit  function  defining  shock  surface 

Entropy 

Transformed  time,  t 
Nondimensional  time 
Time  reference,  t^,  =  r^/-^^ 

Nondimensional  velocity  components  in  9,  (i  space 

Nondimensional  velocity  components  in' the  auxiliary 
F  reference  frame 
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W 

Wn 

X 

a 


y 

At,A?.A0,  A(/) 
5 


n,  ? 

e 

G 

®max 


<}> 


Velocity  components  in  physical  space,  cm/sec 
Shock  velocity  in  radial  direction,  cm /sec 
Shock  velocity  component  in  5  direction 
Dummy  functional  argument 
Angle  of  attack,  deg 

Ratio  of  specific  heats  =  1.4 

Grid  spacing  in  computational  space 

Distance  along  r  between  shock  and  body 

Nondimen sional  transformed  radial 
coordinate  =  (r  -  B)/5 

Vectors  defining  an  auxiliary  cartesian 
coordinate  system 

Transformed  azimuth  angle  =  9 

Coordinate  azimuth  angle,  radians 

Maximum  azimuth  angle  as.  measured  from  wind 
vector 

Flow  density,  gm/cc 

Free- stream  flow  density  on  stagnation 
streamline,  gm/cc 

Transformed  coordinate  meridian  angle  =  4 
Coordinate  meridian  angle,  radians 


SUBSCRIPTS 

1,2 

i,  j,k 
r,  0,  0 


Conditions  fore  and  aft  of  shock  in  Rankine 
Hugeniot  relations 

Free- stream  stagnation  streamline  value 
Indices  for  grid  in  0,  ^ ?  space 
Physical  coordinate  directions 
Free- stream  value 
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